arXiv: 1501.03265vl [cond-mat.stat-mech] 14 Jan 2015 


Quantum transitions driven by one-bond defects in quantum Ising rings 

Massimo Campostrini/ Andrea Pelissetto,^ and Ettore Vicari^ 

^ Dipartimento di Fisica dell’Universitd di Pisa and INFN, Largo Pontecorvo 3, F56127 Pisa, Italy and 
^ Dipartimento di Fisica di “Sapienza” Universitd di Roma and INFN, Sezione di Roma I, 1-00185 Roma, Italy 

(Dated: January 15, 2015) 

We investigate quantum scaling phenomena driven by lower-dimensional defects in quantum Ising- 
like models. We consider quantum Ising rings in the presence of a bond defect. In the ordered phase, 
the system undergoes a quantum transition driven by the bond defect between a magnet phase, in 
which the gap decreases exponentially with increasing size, and a kink phase, in which the gap 
decreases instead with a power of the size. Close to the transition, the system shows a universal 
scaling behavior, which we characterize by computing, either analytically or numerically, scaling 
functions for the gap, the susceptibility, and the two-point correlation function. We discuss the 
implications of these results for the nonequilibrium dynamics in the presence of a slowly-varying 
parallel magnetic field h, when going across the first-order quantum transition at h = 0. 


PACS numbers: 05.30.Rt,05.70.Jk,64.70. qj, 64.60.an 

Quantum phase transitions [ij are phenomena of great 
interest in many different branches of physics. They 
arise in many-body systems in the presence of competing 
ground states. The driving parameters of the transition 
are usually bulk quantities, such as the chemical poten¬ 
tial in particle systems, or external magnetic fields in 
spin systems. In the presence of first-order transitions, 
bulk behavior is particularly sensitive to the boundary 
conditions or to localized defects, hence it is possible to 
induce a quantum critical transition by changing only the 
parameters associated with the defects or the boundaries. 

In this paper, we discuss an example of this type of 
transitions, considering a quantum Ising ring in a trans¬ 
verse magnetic field with a bond defect. In the ordered 
phase, bulk behavior of the low-energy states depends on 
the defect coupling. One may have a magnet phase, in 
which the gap decreases exponentially, i.e., Al ~ e~‘^^ 
with increasing the size L, or a kink phase, in which the 
lowest states are one-kink states and A/, ^ l/L^. Here, 
we analyze the crossover region between these phases, 
showing the emergence of a universal scaling behavior 
controlled by the defect coupling. We also analyze the 
slow nonequilibrium adiabatic dynamics ii across 
this transition. We obtain general time-dependent scal¬ 
ing laws that generalize to the first-order transition case 
those that characterize the Kibble-Zurek (KZ) mecha¬ 
nism at continuous transitions M. 

We consider Ising rings of size L = 2£ 1 in the pres¬ 

ence of a transverse magnetic field and of one bond defect: 

= .fVm - s ^ .f> - < (1) 

i=—i i——l 

where cr|“^ are the Pauli matrices. We set J = 1, and 
assume 5 > 0. Note that periodic (PBC), open (OBC), 
and antiperiodic (ABC) boundary conditions are recov¬ 
ered for C = 1, 0, and —1, respectively. The bond defect 
generally breaks translation invariance, unless ^ = ±1. 


A continuous transition occurs at g = 1, separating a 
disordered (5 > 1) from an ordered {g < 1) phase [l|. 
In the presence of an additional parallel magnetic field h 
coupled to a first-order quantum transition (FOQT) 
occurs at h = 0 for any g < 1, hence we expect the defect 
to be able to change bulk behavior for any g < 1. This 
is the regime we shall consider below. 

Analytic and accurate numerical results can be ob¬ 
tained by exploiting the equivalent quadratic fermionic 
Hamiltonian which is obtained by a Jordan-Wigner 
transformation [1, . We analyze the dependence of low- 

energy properties on the defect parameter () Q . In par¬ 
ticular, we consider the energy differences 

^L,n = En — Eq, Al = A^p, (2) 

where Aq is the energy of the ground state, and En>i 
are the (ordered) energies of the excited levels. The 
magnetization vanishes due to the global Z2 sym¬ 

metry. Thus, we use the two-point correlation function 
G{x,y) = to characterize the magnetic prop¬ 

erties of the ground state. 

For g < 1, we should distinguish a magnet phase 
(C > —1) and a fcmfc phase {( < —1). The lowest states of 
the magnet phase are superpositions of states with oppo¬ 
site nonzero magnetization (±|cri^^|ilA = imp (neglect¬ 
ing local effects at the defect), where [9| toq = (1 — 5^)^/®. 
For a finite chain, tunneling effects between the states 
|-|-) and |—) lift the degeneracy, giving rise to an ex¬ 
ponentially small gap Al [i 3, El- For example, 0 
Al Ri 2(1 — g^)g^ for (^ = 0 (OBC). An analytic cal¬ 
culation gives 

Al-- - we , w = - (1 + C), (3) 

1-5 9 

for C —>■ — I'*’- The large-L two-point function is trivial, 

^ ^ _ G{xi,X2) . ... 

Gr(xi,X2) = -2-^ 1 (4) 
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for xi X 2 , keeping Xi = Xijl fixed (but Xi ^ ± 1 ). 

The low-energy behavior drastically changes for Q < 
— 1 , in which the low-energy states are one-kink states 
(made of a nearest-neighbor pair of antiparallel spins), 
which behave as one-particle states with 0{L~^) mo¬ 
menta [i|. In particular, for C = — 1 (ABC) we have 


Ar = 


l-g 


+ 0{L-^). 


(5) 


The first two excited states are degenerate, thus = 
Al,i = Ai. For ^ < — 1, the ground state and the first 
excited state are superpositions with definite parity of 
the lowest kink | If) and antikink | fi) states. The gap 
scales as ll|, ll 2 | L~^\ we obtain explicitly 


Ar = 


8 C 5 " 


(l_^ 2 )(l _^)2 ^3 




( 6 ) 


On the other hand, A^^^ for n > 2 behaves as L e.g., 


A 


L,2 = 


TT^ 6(1 — C)< 7 ^ 

(T^ + (1 + 0 ( 1 -5)2 


+ O(L-0. (7) 


The two-point function G(x, y) can be perturbatively 
computed for small 5 , obtaining the asymptotic large-L 
behaviors 


G(a:i,X 2 ) = l-|Ai-A 2 | 
G{xi,X2) = 1 — \Xi — A 2 I 


for C = -1, (8) 

I sin(7rAi) - sin(7rA2)| 
- 1 - 


for ( < —1, where Xi = Xiji. We conjecture (and verify 
numerically) that the above formulas can be straightfor¬ 
wardly generalized to the whole ordered phase 5 < 1 by 
simply introducing a multiplicative renormalization, i.e., 
by replacing G with Gr = Gjm^. 

These results suggest that 0 = “1 is a critical point, 
separating the magnet and kink phases. We now show 
that around 0 the system develops a universal scaling 
behavior. We analytically compute (and verify numeri¬ 
cally) the asymptotic behavior of A^^niC), obtaining the 
scaling behavior 


Az.,„(C)~ Ai(Cc)77„(G), (10) 

C. = ^(C-Cc)A, Cc = -1, (11) 
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for L —>■ 00 keeping the scaling variable Cs fixed. The 
scaling functions Di and D 2 are shown in Fig. [T] [l3| . 
The cusp-like behavior at Cs = 0 is the consequence of 
the crossing of the first two excited states at ( = — 1 . 
Generally, the scaling functions are universal apart from 
normalizations of their arguments. In this case, the nor¬ 
malization of Cs is chosen so that the scaling curves for 
different values of 5 are identical. Notice that, once the 
normalization is fixed by using one observable, universal¬ 
ity should hold for any other observable. Numerical re¬ 
sults for the energy differences A^ „ are shown in Fig. [T] 



FIG. 1: (Color online) We show the scaling functions Dn{G)j 
cf. Eq. (1101) . and numerical data for the ratio Ai,n(C)/Ai((^c), 
for n = 1 (bottom) and n = 2 (top) separated by the dotted 
line. Numerical data clearly approach the ^-independent scal¬ 
ing curves Dn{G) (differences are hardly visible). 


they confirm the scaling behavior m- The asymptotic 
large-L behavior is generally approached with corrections 
of order L~^. 

Other observables satisfy analogous scaling relations. 
The two-point function is expected to behave as 


G{xuX2;C)^mlg{X,,X2;Cs), X, = x/e, (12) 


where mo = (1 — 5 ^)^/®. This scaling ansatz can be 
checked by considering the zero-momentum quantities 

X = ^G(0,a;), = ^'^x^G{0,x) (13) 

X X 

(^ is the second-moment length scale). Eq. (fT^ implies 

X/T« m2/^(C.), (14) 

Numerical data confirm them, see Fig. [51 In the lan¬ 
guage of renormalization-group (RG) theory, the defect 
coupling C plays the role of a relevant parameter at the 
magnet-kink transition, with RG dimension = 1. 

An interesting qnestion is whether there is a quantity 
playing the role of order parameter for the magnet-kink 
transition. This is provided by the center-defect correla¬ 
tion b — limL^.oo G(0,.^). Indeed, 6 > 0 for ^ and 
6 = 0 for C < (Cc. Moreover, we observe the scaling be¬ 
havior [ 1 ^ G(0,£) ^ L“^/f,(Cs), with fb{—oo) = 0 and 
fb{oo) = 00 , see Fig. H 

An analogous magnet-to-kink quantum transition can 
be observed in the Ising chain by appropriately tuning a 
magnetic field 5 , coupled to localized at the bound¬ 
aries. Explicitly, we consider (we assume rj > 0) 




5 2 ^^ 


-5 (Ej -crf^). (15) 
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FIG. 2: (Color online) Estimates of the ratio ^/L, supporting 
the scaling ansatz m. Scaling corrections are only visible 
for L < 100. The dashed lines indicate the values of /^(Cs) 
for (s —>■ Too and (s = 0, obtained by matching the scaling 
ansatz with the behaviors in the different phases: /^(oo) = 
l/%/24 from Eq. Q, f^{—oo) ~ 0.098491 from Eq. and 
/{(O) = l/v^dS from Eq. ([8]). The inset shows the crossing 
point of data for different L implied by Eq. HI- 


o 
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FIG. 3: (Color online) Scaling behavior of the center-defect 
correlation G(0,1'). The data of LG{0,£) approach a scaling 
function of The inset shows G{0,i) as a function of C~Cc- 


The analytic computation of the low-energy spectrum 
identifies a particular value of the boundary field, rjc = 
y/1 — g, separating the magnet and kink phases. In the 
magnet phase ij < rjc we have 


Ar. = 


2 gs\/s 2 - 1 ^ 


s-g 

For g > T]c we have instead 


-I- 0(s ), s = 


(16) 




(17) 


with c = 1 for 77 = 77 c and c = 3 for any 77 > 77c. The 
similar nature of the coexisting phases suggests that their 
asymptotic large-L scaling behavior for C Ri ^c and 77 r; 


77 c is the same, i.e., the two transitions belong to the 
same universality class. This is confirmed by analytic 
and numerical computations, although the comparison of 
the results is not straightforward, as the Ising chain (USD 
breaks the Z 2 symmetry, which is instead preserved by 
the Ising ring O- For example, the gap satisfies scaling 
relations analogous to Eq. (HQD. Explicitly, 


Al.„(77) = Az,(77c)£'„(77s), 77 s = ^ ivc - r])L, 

n "" 

with [15[ En{x) — D 2 n-i{x) for a: > 0 and En{x) = 
D 2 n{x) for a: < 0. The reason of the peculiar mapping 
is related to the different behavior under Z 2 of the two 
models. Consider for instance the kink phase. While in 
the Ising ring the lowest states are superpositions of kink 
and antikink configurations, in the case of model (IlSp the 
parity symmetry is broken by the boundary fields, thus 
only kink states are left. Hence, model (flSl) has only half 
of the states of the Ising ring. Moreover, no degener¬ 
acy occurs at 77 = 77 c so that levels must be smooth at 
the transition point, thereby explaining why the mapping 
between the levels must change at the transition point (in 
the ring case, cusps occur at the transition). 

It is worth noting that, at the critical value g = 1, cor¬ 
responding to the order-disorder continuous transition, 
bulk behavior is independent of the boundary conditions 
or of the presence of defects, hence the magnet-to-kink 
transition only occurs for g strictly less than 1. For in¬ 
stance, the gap at g = 1 behaves as A^ ~ for any ( 
or g. Of course, the prefactor depends on the boundary 
conditions; see, e.g., the known results for PBC, OBC, 
and ABC, summarized in Ref. [li¬ 
lt is interesting to reinterpret our results in the equiv¬ 
alent fermionic picture of models m and (ITKI) . In the 
magnet phase, i.e., for (^ > Q and g < gc, respectively, 
the lowest eigenstates are superpositions of Majorana 
fermionic states localized at the boundaries or on the 
defect Bin. In finite systems, their overlap does not 
vanish, giving rise to the splitting A ^ The co¬ 

herence length Iq diverges at the kink-to-magnet transi¬ 
tions as ^ I lns| ^ gc — g and Iq^ ^ ^ in the 
two models, a behavior analogous to that observed at the 
order-disorder transition g —^ 1 “ where ^ I lag|. 

In conclusion, we have shown that quantum transitions 
can be induced by tuning the boundary conditions or by 
changing lower-dimensional defect parameters, when the 
system is at a FOQT. We have explicitly discussed this 
behavior in the case of quantum Ising rings in a trans¬ 
verse field. If g < 1, this model shows a magnet and 
a kink phase, separated by a quantum transition point. 
In its neighborhood, we can define general scaling laws, 
that are analogous to those that hold at continuous tran¬ 
sitions. The same scaling behavior is also observed in 
the XY quantum ring in which one adds additional bond 
couplings [ 1 ^, and in the quantum Ising chain 
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with opposite magnetic fields at the boundaries. The 
universal scaling behavior is essentially the same and is 
uniquely determined by the structure of the low-energy 
behavior in the two phases. We have characterized the 
scaling variable and computed the scaling functions of 
different observables. These scaling behaviors can be 
straightforwardly extended to allow for a nonzero tem¬ 
perature T, by considering a further dependence on the 
scaling variable TL^ = TLp'. Even though we have dis¬ 
cussed the issue in one dimension, we expect the same 
type of behavior in quantum d-dimensional Ising models 
defined in L'^~^ x M boxes with L ^ M, in the presence 
of a (d — l)-dimensional surface of defects or of opposite 
magnetic fields on the boundaries. 

The size dependence of the low-energy spectrum is rele¬ 
vant for the understanding of the nonequilibrium unitary 
dynamics, as it determines the conditions for a nearly 
adiabatic quantum dynamics Significantly differ¬ 

ent behaviors are expected in the magnet and kink phases 
when we add a time-dependent parallel magnetic term, 

Ht = H - h{u) = hou, (19) 

i 

where to is the time scale of the time dependence. Adi¬ 
abatic evolutions across the FOQT (h(0) = 0) require 
very different time scales to- In the magnet phase we 
must have to > while in the kink phase to > L‘^ 

at Cc (»7c) and to > L® for C < Cc [to > L‘^ for V > Vc 
in model (1151) ] . The dynamics in the magnet phase is es¬ 
sentially equivalent to that of a two-level Landau-Zener 
model with an exponentially small gap [iSlIil. At Cc 
for model o and rj > rjc for model di), dynamics be¬ 
comes similar to that at a continuous transition with a 
dynamic exponent z = 2, since there is a tower of ex¬ 
cited states with = 0(L“^). Model ((T|) for C < Cc 
shows again a low-energy dynamics dominated by the 
two lowest states with a gap of order L~^ . The above 
results suggest that nonequilibrium scaling laws, such as 
those describing the KZ mechanism at continuous tran¬ 
sitions (when systems are ramped across a continu¬ 
ous transition at a finite rate), also hold at FOQTs, when 
the time-dependent h crosses the value h = 0. The slow 
nonequilibrium dynamics of the Hamiltonian (|19l) across 
the FOQT h = 0, and in particular the interplay among t, 
to, C (or rj) and L, thus Al, can be described by a scaling 
theory which extends the equilibrium scaling behaviors 
m and m- For this purpose, we use scaling arguments 
similar to those employed to describe the KZ problem . 

Uh = d + z {d is the space dimension, d = 1 in our ex¬ 
ample) is the effective RG dimension of the parallel field 
h (z = 2 and = 3 in the Ising chain (^1. we can define 
an effective length scale associated with h that scales 
as . The dynamical scaling laws for a system 

of size L can be heuristically derived by substituting h 
with the time-dependent h{t/to) into the scaling combi¬ 
nations and This gives t(t/to)^'^^'‘ = 


and Lj^h = T(tAo)^/^'‘ = , 

with T = These considerations lead 

us to conjecture the scaling behavior for the Ising chain 

(4^)) Ri Too /m {xlL, ijr, tIL^, C) , (20) 

(41^42^) ~ ml fg {x,/L, t/r, t/L^, Cs) , (21) 


with given by Eq. (ED. Analogous expressions apply to 
other observables and to the model (fTSl) . Of course, the 
above nonequilibrium scaling theory should be further 
investigated, to get a thorough understanding of KZ-like 
phenomena at FOQTs. 

These issues may also be relevant in the context of 
quantum computing. Adiabatic algorithms rely on suf¬ 
ficiently large gaps during the variation of the model 
parameters b ringing to the ground state of the desired 
Hamiltonian 24-j2|. Thus, FOQTs, at which the gap 
is exponentially small, represent a hard problem 27H29l|. 
As a simple paradigmatic case we may consider the time- 
dependent Hamiltonian m for g < 1 and ^ = 1. Let 
us assume that we want to adiabatically move from the 
ground state with h = —ho at time t = —Ta to the 
ground state with h = ho at t = Ta- This requires an 
exponentially large time scale, i.e. Ta > Our 

results for the C-dependence of the low-energy proper¬ 
ties suggest a way to overcome this hard problem. In¬ 
deed, instead of changing directly h, one could proceed 
as follows. First, one adiabatically changes the system 
varying the bond defect from C, = 1 to C ^ Co which 
corresponds to adding a further single-bond Hamiltonian 
term Hq = In the presence of a finite 

parallel magnetic field h = —ho, the gap is finite. Then, 
one adiabatically changes h, according to Eq. dUD, which 
now requires a time scale Ta > or L® to adiabatically 
go from —ho to ho. Finally, C is increased again to C = !> 
obtaining the ground state of the original problem in a 
total time that scales with a power, and not with the ex¬ 
ponential, of the size. Therefore, by taking advantage of 
particular sensitivity of FOQTs to defects or boundary 
perturbations, one may overcome the problem of an expo¬ 
nentially slow dynamics, which occurs at FOQTs within 
the more standard approaches. 
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